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We compare the superconducting phase-diagram under high magnetic fields (up to J/ = 45 T) 
of Fei_|.j,Seo.4Teo.6 single crystals originally grown by the Bridgman-Stockbarger (BRST) technique, 
which were annealed to display narrow superconducting transitions and the optimal transition tem- 
perature Tc > 14 K, with the diagram for samples of similar stoichiometry grown by the traveling- 
solvent floating-zone technique as well as with the phase-diagram reported for crystals grown by 
a self-flux method. We find that the so-annealed samples tend to display higher ratios Hci/Tc, 
particularly for fields applied along the inter-planar direction, where the upper critical field Hc2{T) 
exhibits a pronounced downward curvature foUowed by saturation at lower temperatures T. This 
last observation is consistent with previous studies indicating that this system is Pauli limited. 
An analysis of our Hc2{T) data using a multiband theory suggests the emergence of the Fulde- 
Ferrel-Larkin-Ovchnikov state at low temperatures. A detailed structural x-ray analysis, reveals 
no impurity phases but an appreciable degree of mosaicity in as-grown BRST single-crystals which 
remains unaffected by the annealing process. Energy-dispersive x-ray analysis showed that the an- 
nealed samples have a more homogeneous stoichiometric distribution of both Fe and Se with virtually 
the same content of interstitial Fe as the non-annealed ones. Thus, we conclude that the excess of 
Fe, in contrast to structural disorder, contributes to decrease the superconducting upper-critical 
fields of this series. Finally, a scaling analysis of the fluctuation conductivity in the superconducting 
critical regime, suggests that the superconducting fluctuations have a two-dimensional character in 
this system. 

PACS numbers: 74.70.Xa, 74.25.Dw, 74.62.Dh, 74.25. fc 



INTRODUCTION 

The Fei+j,(Te,Se) series is a quite unique supercon- 
ducting system: its end member, FeSe is a supercon- 
ductor with a superconducting transition temperature 
Tc ~ 8 K, [l[ which can be increased up to ^ 14 K 
with the partial substitution of Sc for Te, [H, Q or in- 
creased all the way up to 36.7 K by the application of 
hydrostatic pressure. [j| The other end member of this sc- 
ries Fci+j,Te is not superconducting and instead exhibits 
a simuhancous structural and magnetic phase transition 
from tetragonal to monochnic lattice accompanied by an- 
tiferromagnetism (AFM) near Tm ^ 60 to 70 K (Refs. 
0-0). This AFM structure is distinct from those seen 
in the undoped FeAs based compounds, (sl-liot The AFM 
order in Fei^yTe propagates along the diagonal direction 
of the original undistorted Fe square lattice (H], Q , while 
in FcAs the wave- vector of the spin density wave (SDW) 
propagates along the nearest neighbor direction of the 
original Fe square lattice. 

Density functional theory (DFT) calculations indicate 
that the Fermi surfaces (FS) of both FeSe and FeTe are 
similar to those of the other FeAs based compounds, con- 
sequently they should satisfy the nesting condition for the 
wavevector Q^.tt = (tt, tt). [11| This apparent lack of uni- 



versality would seem to question the exclusive role of the 
interband pairing via Qt^^^t magnetic fluctuations in the 
pairing mechanism for all the iron-based superconduc- 
tors. However, the observation of a spin resonance below 
Tc together with an enhancement of the spin fluctuations 
near it, indeed suggests a spin-fluctuations mediated su- 
perconducting pairing mechanism for this system. 12, l3| 
Initial angle resolved photoemission (ARPES) measure- 
ments in the Fei+j,Tc revealed a pair of nearly com- 
pensated electron-hole FS pockets, with no evidence for 



the FS nesting instability associated with Q-^^-n- |14 | 
This could suggest the possibility of a distinct mecha- 
nism for the superconductivity and the magnetic-order 
in the Fe chalcogenides when compared to the Fe ar- 
senides. But a recent ARPES study on superconduct- 
ing Fc1.03Tco.7Seo. 3 reveals a holelike and an electron- 
like FS located at the center and the corner of the Bril- 
louin zone, respectively. These FSs are nearly nested for 
QiT.TT- |15| The same study reports an isotropic super- 



conducting gap along the holelike FS with a gap A of 
^ 4 meV (2A = ksTc ^ 7), thus providing evidence for 
strong-coupling superconductivity. But, this is in con- 
trast with a recent angle-dependent heat capacity study 
in Fei+.ySeo.45Teo.55 which finds evidence for a significant 
gap anisotropy on the electron pockets. [l6j Supercon- 
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ductivity in Fei+^Se was shown to be extremely sensi- 
tive to stoichiometry. [13] In the Fei+y (Te,Se) series, the 
crystal structure resembles that of the iron arsenides [ij 
with the Fe square planar sheets forming from the edge- 
sharing iron chalcogen tetrahedral network. But it also 
allows the partial occupation of iron in the interstitial 
sites of the (Te, Se) layers, resulting in a nonstoichio- 
metric composition for the Fei+j,(Te,Se) series, where y 
represents the excess Fe at the interstitial sites. 0, Q It 
remains unclear how the geometry of the FS evolves with 
the incorporation of interstitial Fe, but it is claimed to 
suppress superconductivity, [l^ 

One remarkable feature of these compounds is their ex- 
tremely large upper critical fields Hc2- For instance, ac- 
cording to Ref. llSl for a Fe1.05Teo.89Seo.11 single-crystal 
with a middle point Tc of just 11 K, one observes an 
Hc2iT ^ K) of - 35 T. Or for a single-crystal of 
Fei.iiTeo.6Seo.4 with a Tc close to its optimum value of 
14 K one obtains iJc2(T ^ K) of - 40 to 45 T. 
Compare these values with those of, for instance, MgB2 
single crystals: = 39 K with H^^{0) = 3.5 T and 
'^c2'(0) = ^"^ T, for fields along the c-axis and along a 
planar direction, respectively. 21[ In fact, to achieve 



in MgB2 upper-critical fields as large as those observed 
in the Fci+y(Te,Se) series, strong impurity scattering is 
introduced to optimize the relative strengths of intra- 
band scattering in a and tt bands of MgB2 . 22 , 2^ The 
significant differences between these two multi-band su- 
perconducting families are attributable to the fact that 
the Fei+y(Te,Se) superconductors are mostly Pauli lim- 
ited [19| while the i?c2S in MgB2 are mostly determined 
by the orbital pair-breaking effect. However, the shape 
of the Hc2 (T) curves for Fei+y(Te,Se) series measured 
in Refs. [20 and [l^ ^-re quite distinct, with the former 
presenting a nearly linear dependence of Hc2{T) for H||c 
and the second displaying a pronounced concave down 
curvature followed by saturation at lower temperatures 
(as expected for the Pauli limiting effect). Such a curva- 
ture has also been reported in Ref. HJ. 

The Pauli-limited behavior of Hc2{T) for Fe- 
chalcogenides results from the fact that they are semi- 
metals with rather low carrier density and low Fermi en- 
ergies -20-50 meV for different FS pockets [3, d, 
unlike the conventional metallic superconductor MgB2. 
As a result, chalcogenides have very short coherence 
lengths ^0 — Hvp/'^ttTc < 1 nm and thus extremely high 
orbitally-hmited H°^^ ~ ^o/Stt^^ ^^ere vp is the Fermi 
velocity in the ab plane, and 0o is the flux quantum. 
This in turn greatly enhances the Pauli effects quanti- 
fied by the Maki parameter um = V^H°2'' / Hp where 
HJT] = 1.84:Tc[K] is the Clogston paramagnetic limit 
[25| . Chalcogenides typically have aM ^ 1 above the 
critical value at which the FFLO instability develops. 
Moreover, because of their very short values of ^0 , chalco- 
genides are naturally in the clean limit £ 3> ^0 which is 
one of the conditions of the FFLO instability, where £ 



is the mean free path due to elastic impurity scattering 
251 . These features of chalcogenides make them good 
candidates to study exotic effects at high magnetic fields, 
in particular, the interplay of orbital and paramagnetic 
pairbreaking for multiband pairing and the FFLO state 
at low temperatures. The small Fermi energies also make 
superconducting properties very sensitive to the doping 
level, allowing one to tune Hc2 by small shifts of the 
chemical potential. 

In contrast to other reports emphasizing on the de- 
tails concerning the synthesis of single crystals [2^, in 
this work we compare the high-field phase-diagram of 
the Fei+y(Te,Se) series, particularly the diagram for the 
optimally doped compound Fei+j,Tei_^~o.6Sea;~o.4 syn- 
thesized by two methods. The first one is based on a 
traveling-solvent floating zone growth technique (TSFZ), 
which leads to crystals of acceptable crystallinity display- 
ing "non-metallic" resistivity, optimal TcS, and transi- 
tion widths ATc ~ 1 to 3 K. We compare their behavior 
with crystals resulting from the Bridgman-Stockbarger 
(BRST) technique which in our case leads to crystals of 
poorer crystallinity (i.e. larger mosaicity), wider tran- 
sitions and frequently to non-optimal TcS in as grown 
crystals, as previously reported by other groups. 
These last single-crystals were subjected to an annealing 
procedure which lead to metallic resistivity, considerably 
sharper resistive transitions ATc — IK, and to supercon- 
ducting transitions comparable to those reported in the 
literature for high quality BRST-grown single crystals as 
measured by SQUID magnetometry. [i^ The annealed 
crystals having a clear metallic behavior preceding super- 
conductivity, display considerably higher upper-critical 
fields when compared to those showing a poor metallicity, 
particularly for magnetic-fields applied along the inter- 
layer direction. This increase leads to a distinct shape of 
the Hc2(T), i.e. from a app roximately linear in T behav- 
ior as reported in Ref. l20l . to the marked concave-down 
curvature followed by saturation at lower temperatures as 
seen in Ref. ^-The analysis of our experimental data us- 
ing a multiband theory of Hc2{T) [2^ shows that Hc2(T) 
is indeed strongly Pauli limited aM — 7 — 10, predicting 
the FFLO state below « 5 K. A detailed single crystal 
x-ray analysis reveals that the annealing process does not 
affect the crystallinity or the degree of mosaicity of our 
single crystals. Nevertheless, our dispersive x-ray analy- 
sis indicates that the annealing process leads to a more 
uniform distribution of interstitial Fe. Annealing yields 
a considerably more isotropic phase diagram which com- 
bined with the enhancement in (^) i suggests that the 
variations in the content of interstitial Fe contributes to 
the suppression of superconductivity. Finally, we show, 
through a scaling analysis of the fluctuation conductivity 
in the neighborhood of the superconducting transition, 
that the superconducting fluctuations in this system are 
two-dimensional in character. 
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SAMPLE PREPARATION 

Single crystals of Fei+xT^Si-ySe.y, 0.05 < a; < 1 and 
0.1 < 2/ < 0.5 were grown by using the traveling- 
solvent floating-zone growth technique (TSFZ) and the 
Bridgman-Stockbarger (BRST) techniques. Starting ma- 
terials with nominal purities of 4N for Fe and 5N for Te 
and Se were used. These were handled in an argon-filled 
glove box which kept the oxygen content below Ippm. 
Mixtures of Fc, To, and Se were ground and pelletized 
and sealed in quartz ampoules under vacuum and heated 
at 400 °C for 24 h. The reacted material was rcground 
and doubly sealed in two quartz ampoules. For samples 
grown through the TSFZ technique, the doubly sealed 
quartz ampoules were loaded into an optical floating-zone 
furnace (NEC Nichiden Machinery 15HD), equipped with 
two 1500 W halogen lamps. The ampoule was rotated at 
20 rpm, and translated at a rate of 1 to 2 mm/h. The 
as-grown samples were annealed by heating to 800 °C for 
48 h, and then cooled at a rate of 100 °C/h to 420 °C and 
held at this temperature for additional 30 h, followed by 
the cooling of the furnace to room temperature. 

A double wall quartz ampoule was also used in the 
case of crystals grown by the BRST technique. The inner 
ampoule was narrowly tipped at its bottom and sealed 
at the top. This Bridgman ampoule was inserting into 
a three-zone temperature gradient vertical furnace and 
lowered at a speed of 4 mm/h. The growth tempera- 
ture was decreased at a rate of 3 °C/h down to 440 °C 
and subsequently quenched to room temperature. Some 
of the crystals as well as Se powder were placed into 
two different quartz crucibles, and both placed inside an 
evacuated quartz ampoule which was sealed under vac- 
uum and heated slowly to 400 °C and then annealed for 
ten days. 

The so-obtained crystals were checked by both single- 
crystal and powder x-ray diffraction (XRD) techniques. 
From as-grown single crystals, shards were cleaved to 
obtain suitable samples of approximately 0.3 mm x 0.3 
mm X 0.03 mm. An Oxford-Diffraction Xcalibur-2 CCD 
diffractometer with Mo Kq, source was used to collect 
reflections. CrysAlis Pro 171.33.55 software (Oxford 
Diffraction) was used for the unit cell refinement and 
analytical absorption correction. Final structure refine- 
ments were conducted using SHELXTL. For the powder 
diffraction runs, the single-crystals were ground by hand 
for approximately 1 minute using an agate mortar and 
pestle. Patterns were detected by a Huber-Guinier 620 
camera with a Rigaku Ultrax Cu direct drive rotat- 
ing anode source. WinPow software was used for unit 
cell parameter least-squares calculations and peak fitting. 
Further details are provided in Ref. lidi 

Scanning electron-microscope energy-dispersive x-ray 
analysis (SEM-EDX) was performed by using a JEOL 
5900 scanning electron microscope (30 kV accelerating 
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FIG. 1. (color online) Top panel: Resistivity as a function of 
temperature for a Fe1.13Seo.45Teo.55 single crystal synthesized 
by the traveling floating zone method under zero magnetic 
field. Inset: Resistive superconducting transition under zero- 
field in a limited temperature range. 



voltage) equipped with IXRF energy dispersion spec- 
troscopy software (IXRF Systems, Inc.) in order to de- 
termine the elemental compositions of the single crystals. 
We found that this EDX set-up would show a tendency 
to underestimate by nearly 3 % the overall fraction of 
Fe when comparing with detailed single-crystal x-ray re- 
finements. We must clarify that the EDX set-up was not 
calibrated against an absolute standard. 

Longitudinal resistivity was measured using a conven- 
tional four-contact method in continuous magnetic fields 
up to 45 T by using either a PPMS or the hybrid mag- 
net of the National High Magnetic Field Laboratory in 
Tallahassee. 

SAMPLE CHARACTERIZATION 

Figure 1 shows a typical trace of resistivity p as a func- 
tion of the temperature for a single crystal grown by the 
traveling-solvent floating-zone method described above, 
whose stoichiometry was extracted by a detailed single 
crystal x-ray analysis refinement. Although the width 
ATc ^ 1 K of the superconducting transition (see inset) 
is relatively narrow when compared with the width of 
the transitions shown for instance in Ref. [l9| . the be- 
havior of the metallic state shows the typical negative 
slope dp/dT < or logarithmic divergence of the resis- 
tivity, which was attributed to the high content of inter- 
stitial Fe and which leads to weak localization, [l^l The 
same characteristic temperature dependent resistive be- 
havior is observed also in as-grown samples synthesized 
by the Bridgman-Stockbarger method (red line), as seen 
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FIG. 2. (color online) Top panel: Resistivity as a function 
of temperature for an as-grown Fei.iiSeo.4Teo.6 single crys- 
tal, synthesized through the Bridgman-Stockbarger method 
(red curve) and under zero magnetic field. After annealing, 
the sample no longer displays the weak-localization like tem- 
perature dependence. Instead it displays a metallic character 
(black trace). Inset: Superconducting resistive transition for 
the annealed crystal under zero field and in a limited tem- 
perature range. Bottom panel: Magnetic susceptibility as a 
function of T for the same annealed single-crystal and respec- 
tively, for magnetic field-cooled (red markers) and zero-field 
cooled conditions (black markers). 



in the top panel of Fig. 2. As seen, the as-grown samples 
tend to exhibit lower superconducting transition temper- 
atures and broader transitions. The annealing procedure 
described above leads to a slightly smaller resistivity at 
room temperature (black line) but most importantly to a 
quite distinct temperature dependence for the resistivity, 
which below T ^ 175 K exhibits a clear metallic depen- 
dence with dp/dT > 0. The magnetic susceptibility on 
the other hand, shows a very clear diamagnetic signal 
(black markers) below the onset of Tc at T = 14 K, as 
seen in the lower panel of Fig. 2. At higher tempera- 
tures, the magnetic susceptibility is virtually tempera- 
ture independent, as reported in Ref. |3l| , and only at 



low temperatures one observes a very mild upturn in the 
susceptibility suggesting the presence of a small amount 
of localized magnetic moments. Although none of the 
susceptibility traces in Ref. |3lJ show a clear diamagnetic 



signal below ~ 14 K but below - 10 K. 
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FIG. 3. (color online) Powder x-ray diffraction pattern for 
two Fe1.11Seo.4Teo. 6 single crystals grown through the Bridg- 
man technique, collected by using the Mo-Kq, line. Black line 
corresponds to the spectrum of an as-grown crystal, while the 
red line is the spectrum obtained from a subsequently an- 
nealed single-crystal. In both cases well-defined Bragg peaks 
are observed and indexed within the space-group Pi/nmm. 
We find no evidence for impurity phases or a change in struc- 
ture/stoichiometry induced by the annealing process. 



We also characterized our single crystals through 
single-crystal and powder x-ray diffraction measure- 
ments. Figure 3 shows two typical powder x-ray spec- 
trum for one as-grown BRST crystal (black line) with 
a nominal stoichiometry Fei.iiSeo.4Teo.6, and for one 
annealed BRST single crystal (red line) from the same 
batch. These "powder" patterns were produced from in- 
tegrating the intensities from single-crystal x-ray diffrac- 
tion measurements using an area detector. This pow- 
der like pattern was obtained on the same single crystals 
used for the magnetic and transport measurements shown 
throughout this manuscript, by superimposing several of 
the 9 to 29 scans along the different crystallographic ori- 
entations. As clearly seen, for both crystals we observe 
the same set of well-defined Bragg peaks, where all peaks 
can be indexed within the tetragonal unit cell having the 
space group symmetry PA/nmm. Most importantly, we 
do not detect the presence of any impurity phase in either 
crystal, concluding that the annealing process preserves 
the original crystallographic structure and therefore the 
stoichiometry of the crystals. From this first analysis it 
would seem that both crystals are equivalent in struc- 
tural quality despite their marked differences in physical 
behavior as seen in the top panel of Fig. 2. 

To clarify if disorder and mosaicity is at the origin of 
these differences in physical behavior, we performed a 
detailed x-ray rocking curve analysis in both single crys- 
tals around specific Bragg peaks. Figure 4 (a) shows the 
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FIG. 4. (color online) (a) X-ray intensity for an as-grown BS 
single- crystal as a function of 29 for an angular scan around 
the Bragg-peak observed along the real-space direction (abc) 
= (200). Red line is a fit to a Gaussian function which yields 
a full-width at half-maximum of ^ 0.5°. (b) X-ray rocking 
curve around the Bragg peak at (200) as a function of the 
angle Q, which corresponds to an arc along a direction nearly 
perpendicular to 9. One observes a very broad feature which 
can be adjusted to three Gaussian distributions (blue lines) 
of crystallite orientations, (c) Same as in (a) but for a Bragg 
peak centered at the (002) position. Red line is a fit to a 
Gaussian having a full-width at half-maximum of ~ 0.25°. 
(d) Same as in (b) but for the peak around the (002) position, 
although in this case two Gaussian distributions are needed 
to adjust the observed broad peak. 



(200) Bragg reflection. By fitting the intensity as a func- 
tion of 29 to a Gaussian distribution, one obtains a full 
wifith at half maximum (FWHM) of ^ 0.5°, suggesting 
at first glance a modest but sizeable mosaic spread among 
crystallites. However, an exploration of the width of this 
peak along an arc 51 whose orientation is nearly perpen- 
dicular to 29 reveals a very broad feature as seen in Fig. 
4 (b). This broad maximum can be adjusted (red line) to 
3 Gaussian distributions (blue lines), with one of them 
showing a FWHM as large as > 2.5°. Each Gaussian 
would represent an ensemble or distribution of crystal- 
lites having a similar relative orientation. The observa- 
tion of several distributions with varying widths indicates 
a large degree of mosaicity around the (200) position, or 
within the planes. A similar study around the (002) po- 
sition shown respectively in Figs. 4 (c) and 4 (d), reveals 
seemingly a smaller degree of mosaicity among crystal- 
lites within distinct planes. 

To clarify if the annealing process affects the structural 
degrees of freedom, for instance, by releasing strain in- 
duced in the single-crystals by the synthesis process, we 
have attempted a similar detailed structural analysis in 
an annealed single crystal of Fei-|-yTeo.6Seo.4. Figure 5 
(a) shows the intensity of the Bragg peak located along 



FIG. 5. (color online) (a) X-ray intensity for an annealed 
BRST single-crystal of Fei-|-yTeo.6Seo.4 as a function of 29 
for an angular scan around the Bragg-peak observed along 
the real-space position [abc] — (200). Red line is a fit to 
a Gaussian function which yields a width at half maximum 
of ~ 0.4°. (b) X-ray rocking curve around the Bragg peak 
at (200) as a function of the angle i}. One observes a very 
broad feature which can be adjusted to at least two Gaussian 
distributions (blue lines) of crystallite orientations, (c) Same 
as in (a) but for a Bragg peak observed along (002) position. 
Red line is a fit to a Gaussian which yields a full width at 
half- maximum of ~ 0.3°. (d) Same as in (b) but for the peak 
around the (002) position. 



the (200) direction as a function of the angle 29, while 
Fig. 5 (b) shows the intensity of the same peak as func- 
tion of the angle il. The first measurement leads to a 
Gaussian peak (red line) having a FWHM of ~ 0.4° while 
the second one can be adjusted to two Gaussian distri- 
butions (blue lines) with the widest distribution having 
a FWHM of ~ 1.6°. Fig. 5 (c) shows the intensity of the 
Bragg peak observed as the angle 29 is scanned around 
the (002) direction. One extracts a FWHM of ~ 0.3° 
by fitting its intensity to a Gaussian distribution (red 
line). Finally, Fig. 5 (d) shows the intensity of the Bragg 
peak observed at the (002) position as a function of the 
angle 51. As shown in the figure, to reproduce the ob- 
served broad feature one would need to fit it to at least 
three Gaussian distributions with the widest one having 
a FWHM of > 2.1°. 

From this detailed comparison between the as-grown 
and the annealed BRST Fei-|_yTeo.6Seo.4 single crystals, 
we can firmly conclude that the annealing process has 
no or perhaps just a marginal effect on the degree of 
crystallinity of the BRST grown samples. Although as 
seen in Fig. 2 (a), the behavior of the metallic state, 
and even the width of the superconducting transition is 
affected by this treatment and as we shall see below, 
it also affects Hc2{T). In Figs. 6 (a) through (d) we 
show a similar rocking curve analysis for a typical TSFZ 
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FIG. 6. (color online) (a) X-ray intensity for a single-crystal of 
Fei+j;Teo.55Seo.45 grown by the traveling floating zone method 
as a function of 29 for an angular scan around the Bragg- 
peak observed along the (200) position. Red line is a fit to a 
Gaussian function which yields a width at half maximum of 
--^ 0.2°. (b) X-ray rocking curve around the Bragg peak at 
(200) as a function of the angle Q. One observes a very broad 
feature which can be adjusted to a superposition of least ten 
Gaussian distributions (blue lines) of crystallite orientations, 
(b) Same as in (a) but for a Bragg peak observed along (002) 
position. Red line is a fit to a Gaussian which yields a full 
width at half- maximum of ~ 0.3°. (d) Same as in (b) but for 
the peak around the (002) position. 



y) throughout the sample. 

However, this observation does not address the ori- 
gin of for instance, the "non-metaUic" hke resistivity ob- 
served in the TSFZ-grown crystals when compared to the 
metallic response seen in annealed BRST-grown crystals. 
For the first ones, we were able to observe clear Bragg 
peaks and proceed with single-crystal x-ray refinements 
in shards taken from them, indicating relatively large val- 
ues for y ranging from 0.07 to 0.12. The crystallographic 
information file for one of our single-crystals is avail- 
able as supplemental material [s^. We also performed 
EDX analysis on the TSFZ single crystals, finding simi- 
lar amounts of excess interstitial Fe, and typical values for 
Ay ~ 2 %. Remarkably, for the second ones and despite 
the physical characterization shown above suggesting rel- 
atively high-quality single-crystals, we were unable to ex- 
tract well-defined integrated intensities to proceed with 
single-crystal x-ray refinements. For these crystals, EDX 
indicates nearly stoichiometric values for y, i.e. y ^ 
although, as previously mentioned, it tends to underesti- 
mate the content of Fe when compared to single-crystal 
x-ray refinements. These observations indicate that the 
physical properties of the 11 Fe-chalcogenides are nearly 
oblivious to their degree of mosaicity but depend on the 
amount of interstitial Fe and on how evenly distributed 
it is throughout the sample. 

SUPERCONDUCTING PHASE DIAGRAM AT 
HIGH FIELDS 



Fe1-1-yTeo.55Seo.45 single-crystal, and as can clearly be 
seen in Figs. 6 (b) and (c) the degree of mosaicity in 
this single crystal, is higher (spans over 6°) than what 
we found for the BRST crystals. Nevertheless, the width 
of the superconducting transition (sec Fig. 1) is consider- 
ably sharper than the width of the transition in as grown 
BRST crystals (see Fig. 2(a)) Therefore, from our study 
we can state that the mosaicity, or the relative orientation 
between stacked planes in a single-crystal, has virtually 
no effect on the superconducting properties (such as the 
width of the transition) of the 11 series. 

Nevertheless, by comparing the EDX data taken on 
several annealed and non-annealed BRST-grown single- 
crystals, on approximately a dozen of collection sites in 
each crystal, one notices marked differences between an- 
nealed and non-annealed crystals: in average the non- 
annealed crystals displayed a variability Ay of ~ 7.7 % 
in the fraction y of interstitial Fe and a Ax ~ 8.8 % 
in the fraction x of Se, when compared respectively to 
Ay = 1.8 % and Ax = 6.1 % for the annealed samples 
(Ay is the average of the fiuctuations with respect to the 
average value of y). EDX indicates that the annealing 
process has virtually no effect on the average value of ei- 
ther X or y, so it just homogenizes both values (specially 



In this section we compare the high-field supercon- 
ducting phase-diagrams of several crystals synthesized 
by either the floating zone or the Bridgman-Stockbarger 
technique. In each case, only the crystals display- 
ing the sharpest superconducting transitions ATc = 
T(p(0.9rj - T{p{0.1Tc)) < 1.0 K (where pi0.9Tc) and 
p(O.lTc) are the values of the resistivity at respectively, 
the onset of the resistive transition or 90 % of p, and the 
foot or 10 % of p at the transition) were chosen for the 
electrical transport measurements under high fields. As 
indicated by Fig. 2, for the crystals grown by BRST the 
method, only the annealed samples showed the optimum 
Tc and the sharpest superconducting transitions among 
all of the measured samples. 

Figures 7 (a) and (b) show the resistance as a func- 
tion of temperature for an annealed Fei-|-j,Teo.6Seo.4 sin- 
gle crystal for different fields H||a6 and H||c, respectively. 
Clearly, fields along the c-axis produce broader transi- 
tions as a function of temperature. Figures 7 (c) and (d) 
show the resistance as a function of magnetic fields H||a6 
and H||c (up to 77 = 45 T) for the same single crystal for 
different temperatures. The transitions become progres- 
sively sharper in field as the temperature decreases, con- 
sistent with the reduction of vortex fluctuations at lower 
T. We have also collected similar data-sets for several 
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FIG. 7. (color online) Left panel: Resistance as a function of 
temperature T across the superconducting transition for an 
annealed Fei.iiSeo.4Teo.6 single crystal (nominal stoichiom- 
etry) synthesized by the Bridgeman technique, and respec- 
tively for fields along an in-plane direction (a) and the inter- 
plane c-axis (b). Right panel: Resistance as a function of 
field H for the same single-crystal and for different temper- 
atures, and respectively for an in-plane direction (c) and the 
inter-plane c-axis (d). 



single-crystals grown by the TSFZ-method with differ- 
ent Se and Fe contents but the corresponding raw data is 
not shown here. The stoichiometric fractions of Fe and Se 
were determined through a detailed single-crystal x-ray 
diffraction refinement of small shards cleaved off from 
each single-crystal. For these TSFZ samples the mea- 
sured superconducting transition temperatures ranged 
from Tc^l2K{y = 0.26, x or fraction of Se = 0.16) to 
the optimal Tc value > 14 K (y = 0.06 and x = 0.44). 

In Fig. 8 (a) we compare the superconducting 
phase-diagram at high-fields for our annealed BRST 
Fei+j,TeQfiSeo.4 single-crystal (circles) with the results 
of Ref. |20| (squares), and in Fig. 8 (b) with the phase- 
diagram obtained for several TSFZ single-crystals (var- 
ious markers excluding circles). Here, all points corre- 
spond to the middle point of the resistive transition. In 
either graph, both the field- and the temperature-axis 
are re-scaled with the zero-field of each sample (cor- 
responding to the middle point of the zero-field resistive 
transition). Here, we use this scaling for comparative 
purposes only, and not to indicate that these compounds 
are placed within the dirty limit. The plot in Fig. 7 (a) 
leads to a similar scaling quality when compared with 
a Hc2{T)/T'^ as a function of T/Tc plot, previously ob- 
served in some of the Fe arsenide compounds |35[ and 
which was deemed to be consistent with an effective clean 
limit scenario due to the small superconducting coher- 
ence lengths of the Fe-pnictide compounds. As clearly 
seen in both figures, the annealed BRST-grown crystal 
has considerably larger upper-critical fields than all the 



FIG. 8. (color online) (a) Comparative superconducting 
phase-diagram under high magnetic fields where the upper 
critical field Hc2 normalized by the superconducting transi- 
tion temperature Tc, is plotted as a function of the reduced 
temperature t = T /Tc for two single crystals: respectively for 
the Fei+aSeo.4Teo.6 single crystal in Ref [53 and for our an- 
nealed Fe1.11Seo.4Teo. 6 (nominal stoichiometry) BRST single- 
crystal, (b) Superconducting phase-diagram for several TSFZ 
grown single- crystals, where the phase-diagram for our an- 
nealed Fe1.11Seo.4Teo. 6 BRST single-crystal is also included. 
Open and solid markers depict the phase-boundary between 
metallic and superconducting states for fields along the c-axis 
and along a planar direction, respectively. All points in both 
graphs were extracted from the middle point of the resistive 
transition. 



other crystals, particularly for fields along the c-axis and 
at intermediary reduced temperatures, i.e. from ~ 35 to 
- 60 % larger iJc2S for T/T^ ~ 0.5. For T/T^ ~ 0.2, the 
upper-critical field still is ~ 5 % larger than the value 
reported in Ref. [2^ On the other hand, the anisotropy 
parameter 7 ~ 1.8 for T Tc is slightly smaller than the 
values reported in Ref. (7^ 2 for T -> T^), Ref. [H 
(7 ~ 3 for T ^ Tc), and Ref. [H (7 ~ 3.5-4 for T TJ. 
Given that the level of structural disorder (i.e. mosaicity) 
seen in crystals showing quite broad transitions and lower 
TcS (as-grown BRST crystals) is very similar to the level 
of mosaicity seen in the annealed crystals which display 
relatively sharp transitions, optimum Tc, and higher i?c2, 
one is lead to conclude that the superconducting phase 
diagram of the Fei+yTei_2,Sea; series for a given value 
of X is mostly controlled by Ay and Ax respectively the 
overall fluctuation on the value of y the fraction of inter- 
stitial Fe and in the value of x the fraction of Se. 

Another intriguing observation is that the supercon- 
ducting phase-diagram for the more Fe stoichiometric 
samples shows a pronounced downward curvature as T 
is lowered which wa s attrib uted by several authors to the 
Pauli limiting effect 



19ll33|. in contrast with for example. 



the linear dependence seen in Ref. |20D for H^2- Indeed, 
all Hc2/Tc ~ 3 ratios shown in Fig. 7 are above the BCS 
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mogencities. The fit yields Hf = (72.3 ± 3.5) T and 
hI"^ = (62.6 ±2.4) T for the Pauli limiting field, respec- 
tively for fields along the c-axis and along the ffl6-plane, 
and hI'' = (160. 1± 17.3) T and H^^ = (439±136) T for 
the orbital hmiting fields wliicli are in very good agree- 
ment witli those reported in Ref. [3^. Defining tlie effec- 
tive GL coherence lengths, ^ab{T) = ^afc(l - T/TJ^^/^ 
and ^c(r) = Cc(l - T/Tc)"^/^, we obtain a rather small 
in-plane value of ^ab ^ {(jJo/'^TrHty/'^ = 14 ± 4 A and 
an even smaller c-axis value of — S,abHo'" / Ho°'^ ^ 5.1 
A, which is shorter than the inter-planar distance c ^ 6 
A 34 1, in agreement with Ref. [ssl . 



FIG. 9. (color online) Upper critical fields Hc2 as a function 
of 1 — t where t — T/Tc is the reduced temperature for an an- 
nealed Bridgman-Stockbarger grown single-crystal. Red lines 
are fits to Eq. 1. It yields values for the coherence lengths ^ab 
and which are in good agreement with those of Ref. Issl . 



paramagnetic limit, Hc2/Tc = 1.84, indicating that the 
Pauli pairbreaking is very essential. 



Two-band analysis of the Hc2 data 

In order to evaluate the contributions of both orbital 
and Pauli pairbreaking effects for either field orientation, 
and to compare the results with those extracted from 
heat capacity measurements, we first analyze our Hc2{T) 
data at temperatures close to Tc where the Ginzburg- 
Landau theory yields Ref. [2^, [s^ : 



H 



H _^ T 
Ho Tc 



(1) 



Very close to the critical temperature, (Tc — T)/Tc <C 
(Hp/Ho)^, the first paramagnetic term in the left hand 
side is negligible and Eq. ([!]) yields the orbital linear GL 
temperature dependence, Hc2 = Ho{l — T/Tc). At lower 
temperatures, (T^ - T)/Tc > {Hp/Haf, the Pauh limit- 
ing field Hp dominates the shape of Hc2iT) oc (1 — t)^/^ 
even in the GL domain if Hp < Hq. The latter in- 
equality is equivalent to the condition that the Maki 
parameter um ^ Ho/ Hp > 1 is large enough, assur- 
ing that the paramagnetic effects are essential. Shown 
in Fig. 9 are the log- log plot of our Hc2{T) as a func- 
tion of 1 — T/Tc where the red lines are fits to Eq. 
([1]). These fits are excellent for the high- field region 
(excluding the highest fields), but less so very close to 
Tc, probably due to the relatively large errors in deter- 
mining temperatures (AT ~ 25 mK) which are inherent 
to transport measurements, or perhaps due to broad- 
ening of the resistive transition due to local Tc inho- 
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FIG. 10. (color online) High-field superconducting phase- 
diagram or upper critical-field normalized by the super- 
conducting transition temperature at zero field as function 
of the reduced temperature, for an annealed BRST-grown 
Fei+5Seo.4Teo.6 single crystal, where the red line is a fit to the 
model described in the text. Blue line depicts the tempera- 
ture dependence of the modulation wave-vector for a FFLO 
state, which would explain the shape of the phase-boundary 
at lower temperatures. 



The phenomenological Eq. ([T]) only describes the range 
T fa Tc, but does not take into account the possibility of 
the FFLO state which occurs if um > etc is large enough 
(«£ ~ 1-8 for a single parabolic band |25|). To analyze 
our Hc2{T) data at all temperatures and to reveal the 
microscopic meaning of the scaling fields Hp and Ho, we 
use a two-band generalization of the WHH theory [2^ 
in the clean limit, taking into account both orbital and 
paramagnetic pair-breaking, and the possibility of the 
FFLO with the wave vector Q{T,H). In this case the 
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equation for Hc2 is given by 



aiGi + 026*2 + G1G2 = 0, 



(2) 



Gi =lnt + 2e« 



: tan 



Re^:/ , 

/ u\/b 



due " X 



n + 1/2 V& \ t(n + 1/2) 

Int + 2e'*9'Rc^ / due~''\ 

n -J a\, 

u t 



Go 



n+1/2 



tan^ 



,j f uy/brj 



t{n+ 1/2) + iab 



(3) 



The FFLO wave vector Q{T,H) is determined self- 
consistently by the condition that the solution Hc2{T, Q) 
of Eq. ^ is maximum, ai ~ (Aq + X-)/2w, 02 ~ 
(Ao - X^)/2w, A_ - An - A22, Aq = (A^ + 4Ai2A2i)i/2, 
w = A11A22 - A12A21, t = T/Tc, and 



4:fi(j)ogikBTc TTkBTcTn 



(4) 

s = e2/ei. (5) 



Here i"; is the in-plane Fermi velocity in band / = 1,2, 
ei = mf''/m'i is the mass anisotropy ratio, /i is the mag- 
netic moment of a quasiparticle. An and A22 are the in- 
traband pairing constants, and A12 and A21 are the in- 
terband pairing constants, and a ss 0.56aA/. The factors 
51 = 1 + An + IA12I and 52 = 1 + A22 + IA21I describe 
the strong coupling Eliashberg corrections. For the sake 
of simplicity, we consider here the case of ei = £2 = e for 
which s = 1, and Hat is defined by Eqs. ([2]) and (3) 
with gi = (72 and rescaled q — qe~^^'^, a — >■ ae~^^^ and 
Vb ^ ei/4y^ in and ^ e^/^Vrjb in G2 

Equation ([2]), which describes Hc2{T) at all tempera- 
tures, simplifies close to Tc where it reduces to Eq. ([T]) 
with 



487r. 



■ V. 



2 -r {vi 

2TikBTc 



(6) 
(7) 



where C(3) « 1.202. The slope H'^^ = \dHc2/dT\ is max- 
imum at Tc where Eqs. (2) and ([6]) yield 



Hc2 - 



247r0ofc|rc 



(8) 



For identical bands, {vi = V2 ~ v, An = A22, c+ = C- = 
1/2), Eq. ([U reduces to the single-band GL expression, 
H'c2 = 4>Ql2T:Tcil, where = (S^;/^fcBTe)[7C(3V48]i/2 
is the GL coherence length in the clean limit [29|. For 
the pairing (c+ — >■ c_ — > 1/2), the dependence of H'^j 



on the materials parameters resembles the behavior of 
Hc2(T) in the s"*""*" dirty limit: for strong band asymme- 
try (?7 ^ 1 or ?7 1), Hc2 in Eq. ([8]) is limited by the 
band with larger Fermi velocity, similar to Hc2 mostly 
limited by the band with larger difFusivity for the s"*""*" 
case [11] . Paramagnetic effects decrease the slope of 
and reduce the effect of band asymmetry. 

Figure 10 shows the fit of the measured Hc2{T) to 
Eq. ([2|) for H\\ab. For this field orientation, the resis- 
tive transitions in Fig. 7 (c) are considerably sharper 
than for H\\c. This allowed us to clearly define the mid- 
dle point of the resistive transition for the traces taken 
at lower temperatures, increasing the number of the Hc2 
data points. For the sake of simplicity, we consider here 
the case of ei = £2 = e for which 7 = e^^^^, and Hat is 
defined by Eqs. ([2]) and (3) with rescaled q qe~ 
ae-y2 and Vb -> e^/'^Vb in Gi and y/7jb e^l^^ 



a 



-3/4^ 

in G2 [29|. The fit in Fig. 10 was done for pairing 
with An = A22 = 0, A12A21 = 0.25, 77 = 0.9, H\\ab, 
a = 2.3. Equation Q describes H^2{T), H^i{T) and 
7ff(T) = 6||(T)/Ve5_L(r) where &||(r) is the solutions of 
Eq. © for i7||c very well. 

To see if the observed values of 11^2 are at least qualita- 
tively consistent with Eq. ([8]) and the materials param- 
eters of chalcogenides, we use the ARPES data 
which give Ep 20 - 50 meV. Taking Ep = 30 meV 
in Eq. dH]) yields H'^f [T/K] ~ 0.647(m/TOo), so the ob- 
served H'cf = Hf /Tc = 31.4 T/K and 7 = Hf/H^ « 
2.75 corresponds to the value of m ~ 15too consistent 
with the ARPES data of m/mo = 3-20 for different FS 
pockets 43|. These estimates show that the paramag- 
netic effects for Fe-11 chalcogenides are indeed essential, 
and the parameter defined by Eq. (jH) is greater 

than 1. For Ep = 30 meV and = 14 K, we ob- 
tain a"'' ~ 0.137m/mo, giving a"'' ~ 5. The values of 
a > 1 indicate that paramagnetic pair-breaking becomes 
so strong that it can trigger the FFLO instability at lower 
temperatures. 

As follows from the results shown in Fig. 10, the 
anisotropy parameter 7^f(T) decreases as T decreases. 
This behavior reflects the significant role of the Zeeman 
pair-breaking in FeSeTe given that ~ ce/V^- = 2.3 
for HJ-c is well above the single-band FFLO instabil- 
ity threshold, a « 1 [2^. In this case 7_f/(T) near Tc 
is determined by the orbital pairbreaking and the mass 
anisotropy e, but as T decreases, the contribution of the 
isotropic Zeeman pairbreaking increases, resulting in the 
decrease of jh{T). 

The fit to our experimental data based on the solution 
of Eq. ^ predicts the FFLO transition for H\\ab. In 
this case Fig. 10 shows that the FFLO wave vector Q{T) 
appears spontaneously at T < Tp « 0.35rc ~ 5K, corre- 
sponding in optimally doped samples to magnetic fields 
i7 > 42 T. Here Q{T) increases from zero at T = Tf to 
the maximum Qo at T = 0, where Qo = 47rTc6(0)^/^/ft?;i. 
In addition to the strong Pauli pairbreaking a > 1, the 
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FFLO transition requires weak impurity scattering in the 
clean limit £ 2> where I is the mean free path due to 
elastic scattering on impurities [2^ . The latter condition 
is likely satisfied in the crystals studied in this work, given 
the very short coherence lengths ^ab ~ 14 A and ~ 5 A 
extracted from the fit described above. Thus, the chalco- 
genides are very good candidates to study the FFLO 
transition by magneto-transport, specific heat, magnetic 
torque or NMR measurements. 

ANALYSIS OF THE CRITICAL REGION 

Recent heat-capacity measurements on FcSco.sTco.s 
under field [36[ found that it strongly resembles behav- 
ior previously observed in the much more anisotropic 
high-Tc cuprates. where strong fluctuation effects have 
been found to wipe out the phase transition at Hc2- In 
FeSeo.sTeo.s these fluctuations were found to be strongly 
anisotropic [s^. Here, we proceed to study the criti- 
cal regime preceding the superconducting transition in 
our annealed BRST-grown single-crystals. In the crit- 
ical region of the metal to superconducting transition, 
and in the limit of strong magnetic fields, one expects a 
scaling form for the thermodynamic functions. At very 
high fields if the quasiparticlcs arc confined within the 
lowest Landau level, the transport of charge carriers be- 
comes nearly one-dimensional along the direction of the 
applied field. Fluctuation effects close to the supercon- 
ducting transition are expected to be enhanced by the 
lower effective dimensionality of the system. In the criti- 
cal regime the fluctuation conductivity was calculated in 
Ref. |37i by including a quartic term in the free energy 
within the Hartrce approximation, and obtained a scal- 
ing law for the fluctuation conductivity Act in magnetic 
fields, in terms of unspecified scaling functions F2D and 
F^Y) , valid for two-dimensional and three-dimensional su- 
perconductors, [s^l respectively: 



At(H)2D = ( I 



1/2 



F2D a 



T - T,{H) 



1/3 



(9) 



(10) 



Where a and /3 are characteristic constants of a specific 
material. Presumably, these functionals remain valid 
even if the quasiparticlcs are not confined to just the 
lowest Landau level, but are placed in a few higher Lan- 
dau levels, as long as the interaction between quasi- 
particles in distinct Landau levels remains negligible. 
[H, [3^ In the cuprates, either type of scaling has been 
observed in materials having distinct degrees of elec- 
tronic/superconducting anisotropics, [i^l 

We define the contribution of the fiuctuations to the 
conductivity Act as the difference between the normal 



state conductivity (t„ = l/p„ (p„ is the normal state 
resistivity) and the measured conductivity a = 1/ p. Here 
(T„ is obtained from a polynomial fit of the resistivity 
in a temperature interval AT 2Tc above Tc, where 
the contribution of the fluctuations to a should become 
negligible. Here, Tc is defined as the middle point of the 
resistive transition, or the 50 % 

In the left panel of Fig. 11, we plot Acr(i//r)i/2 ver- 
sus [T - TcXH)]I{THY/^ for the annealed BRST single- 
crystal for both field orientations, as expected for the 
two-dimensional scaling relation in Eq. For fields 

above a relatively small value of just 4 T, and for both ori- 
entations, this two-dimensional scaling clearly succeeds 
in collapsing all the curves over an extended temperature 
range particularly above Tc- We checked that this state- 
ment is still valid if the onset of the resistive transition 
(90 % of the resistance of the normal state) is used as the 
criteria to define T^- Given that the high field resistive 
transitions shown in Fig. 7 are remarkably broad, it is 
not possible to unambiguously define the functional form 
of Pn{H) in order to also collapse those curves in Fig. 
11. To compare this two-dimensional scaling relation 
with the three-dimensional one, in the right panel of Fig. 
11, we plot AaiJi/VT^/s versus [T - Tc{H)]/ {THf/^ 
as expected for the three-dimensional case described by 
Eq. (fTO|) . Clearly, the two-dimensional relationship pro- 
vides a better scaling of the data both above and below 
Tc- Indeed, in the three-dimensional case, the curves at 
different fields open like a fan above Tc and also tend 
to separate from each other below this critical temper- 
ature. On the other hand, in the 2D case, the curves 
above 4 T basically overlap over the entire range. The 
2D behavior of fluctuation conductivity, already observed 
in the Fe arsenides, [4l|, IH] implies that the inter-plane 
coherence length is smaller or in the order of the 
inter-plane spacing as found above, or, equivalently, that 
the coupling between adjacent superconducting planes is 
smaller than the intra-plane condensation energy. Notice, 
that our scaling analysis is consistent with the report of 
very anisotropic superconducting fluctuations according 
to heat capacity measurements [36| . Two-dimensional 
scaling also suggests that such fluctuations contribute to 
the broad transitions seen in Fig. 7, which become con- 
siderably sharper at lower temperatures. 



DISCUSSION 

There are several puzzling aspects in the 
Fei-|-yTei_2;Sea; family that have yet to be understood 
and/or conciliated. For instance, optical spectroscopy 
measurements in Fei.oeTeo. 8380.14 indicates the absence 
of a Drude peak in optical conductivity, suggesting 
the absence of well-defined coherent quasi-particles in 
this material 0. Although, for FeSeo.42Teo.58 the 
in-plane optical conductivity is found to be indeed 
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FIG. 11. (color online) Left panel: Aa{H/T)'^^^ versus AT{H) / {THf^ for the annealed BRST single-crystal and for both 
field orientations, as expected for the two-dimensional scaling relation of Eq. (2). Here, Act = (1/p — l/p„) where 1/p is 
the resistivity across the superconducting transition and 1/pn is the resistivity in the normal state preceding the transition 
(see main text), while Ar(7/) = [T - T^{H)]. Right Panel: AaH^^^ /T'^'^ versus [T - Tc{H)]/ {THf/^ as expected for the 
three-dimensional scaling relation in Eq. (3). 



describable within a Drudc-Lorentz model [47| but with 
an incoherent response along the inter-planar direction 
[48l |. Therefore, one could even ask if the concept of 
Fermi surface would be applicable to these compounds. 
In sharp contrast, a recent angle resolved photoemission 
study on FeSeo.42Teo.58 [S] reveals a Fermi surface 
which is in relative good qualitative agreement with 
the generic Fermi surface proposed for the family of Fe 
pnictides. Namely, composed by two concentric hole-like 
cylindrical surfaces at the F point and two electron-like 
cylinders at the M point of the Brillouin zone, which arc 
characterized by a large, FS-sheet dependent, effective- 
mass enhancement. This indicates that correlations 
are particularly relevant for these materials which is 
consistent with the size of the anomaly observed in the 
heat capacity at the superconducting transition [H^l 
or the large Hc2 reported here and by other groups 
[20I 24 . [sif . This geometry would be consistent with 



the observation of a resonance mode characterized by a 
wave-vector connecting both types of cylindrical Fermi 
surfaces 1^ pointing towards an itinerant, nesting-like 
mechanism for the origin of the spin-fluctuations in 
these systems. Although, this scenario seems difficult to 
conciliate with the rather large magnetic moment of 1.6 
to 1.8 yUs for the non-interstitial Fe(I) ion as extracted 
from careful magnetic susceptibility measurements in 
superconducting samples [Hlj, and does not explain the 
observation of short-range magnetic-order in (supercon- 
ducting samples) for the same wave-vector as the one 



extracted for the magnetically ordered state of Fci+^Te 
0. The interstitial Fe(II) is presumably characterized 
by an even larger moment of ~ 2.5 jiB @, Uli and 
from the superconducting perspective it should act as a 
magnetic impurity. 

A recent DFT + DMFT study suggests that a com- 
bination of Hund's rule and structural degrees of free- 
dom, such as the pnictogen height and the bond angle 
between Fe and the pnictide element, are the key pa- 
rameters defining the role of correlations, i.e. the degree 
of localization and therefore the concomitant size of the 



magnetic moment of Fe in iron pnictide compounds [52 1 . 
For instance, the larger the size of the pnictogen atom or 
equivalently the larger the distance between Fe and this 
element, or the greater the deviation with respect to the 
ideal tetrahedral angle of 109.5°, the larger the degree of 
localization in the Fe site. Or equivalently, the stronger 
the correlations on those Fermi surface sheets having a 
marked t2g character which would favor the development 
of a magnetic instability. It is therefore easy to speculate 
on the role played by the interstitial Fe(II): it should favor 
local lattice distortions that are likely to further distort 
the Fe-pnictogen angle away from 109.5° favoring local 
magnetism in detriment of superconductivity. The larger 
the content in Fe(II) the stronger this detrimental ef- 
fect, while larger fiuctuations in y could create randomly 
distributed nearly magnetic patches having a stronger 
Fe(II) content which could lead to pair-breaking effects. 
This rather simple scenario would explain why the an- 
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nealed BRST-grown samples display clearer and sharper 
superconducting transitions when compared to the non- 
annealed ones despite a similarly poor crystallinity. 

As for the presence of the FFLO state [2^, it is pre- 
sumably sensitive to impurities being difficult to concil- 
iate with the presence of magnetic Fe(II) atoms. How- 
ever, our results suggest that our BRST-grown samples 
(annealed as well as non-annealed) are closer to having a 
nearly stoichiometric Fe content. Furthermore, the low- 
T upturn in H^2 (T) indicative of the FFLO instability as 
observed in this work, was also observed by other groups, 
see for example, Ref. [IJ. Finally, the interplay between 
magnetism and the FFLO-state is still poorly understood 
53|, so its existence in the Fei-f.yTei-xSex will only be 
clarified by future calorimetric, magnetic torque or NMR 
measurements at high fields. 

Finally, given that the inter-plane coherence length S,c 
is in the order of the inter-planar distance it is relatively 
easy to understand why the mosaic spread associated 
with relative orientation among Fe[Se,Te] planes has little 
effect on the superconducting properties of these mate- 
rials. Although the superconducting anisotropy is small, 
the size of implies that the Cooper pairs are virtually 
localized within the planes thus remaining almost oblivi- 
ous to the potential pair-breaking effects imposed by the 
grain boundaries located between the planes. 



CONCLUSION 

In conclusion, we find that the Fe[Te,Se] planes are 
well stacked along the c-axis, forming domains with well 
aligned planes spanning a certain coherence length along 
the c-axis. These domains display a small twist to each 
other, in the order of a degree or less relative, defin- 
ing grain boundaries. Our study strongly indicates that 
these domain boundaries have little or no effect on the 
superconducting properties of the Fei+yTei_2;Sea; series. 
However, a relatively weak disorder in the distribution of 
interstitial Fe is found to be detrimental to the supercon- 
ducting properties of the Fei+yTei-j^Se^; series, broad- 
ening the width of the superconducting transition and 
suppressing the diamagnetic response. We attribute this 
to local lattice distortions associated with the intersti- 
tial Fe(H) which are likely to locally increase the elec- 
tronic correlations leading to a higher degree of carrier 
localization favoring magnetism in detriment of super- 
conductivity. In samples displaying the sharpest super- 
conducting transitions and a very clear diamagnetic sig- 
nal, we observe a mild upturn in the upper critical field 
at lower temperatures consistent with the Fulde-Ferrel- 
Larkin-Ovchinnikov instability. These samples display 
metallic resistivity suggesting also the formation of co- 
herent quasi-particles at lower temperatures in contrast 
to previous reports (46j . and implying that a relatively 
mild excess of randomly distributed Fc(II) atoms may be 



enough to suppress the phase coherence. Finally, our es- 
timates of the inter-plane coherence length ^ 5.1 < 6 
A where c ~ 6 A is the inter- planar lattice constant, is 
consistent with the quasi-two-dimensional superconduct- 
ing fluctuations in the critical region (also observed by 
other groups [4l|, |42|). 
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